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Abstract 

In a recent work, a dual formulation of group field theories as non-commutative quantum 
' field theories has been proposed, providing an exact duality between spin foam models and 

, non-commutative simplicial path integrals for constrained BF theories. In light of this new 

framework, we define a model for 4d gravity which includes the Immirzi parameter 7. It 
reproduces the Barrett-Crane amplitudes when 7 = 00, but differs from existing models 
^ ' otherwise; in particular it does not require any rationality condition for 7. We formulate the 

On| , amplitudes both as BF simplicial path integrals with explicit non-commutative B variables, 

' and in spin foam form in terms of Wigner 15j-symbols. Finally, we briefly discuss the 

correlation between neighboring simplices, often argued to be a problematic feature, for 
example, in the Barrett-Crane model. 
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■ 1 Introduction 

rS 
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■ Group field theories [1] are quantum field theories showing up as a higher dimensional general- 
ization of matrix models in background independent approaches to quantum gravity [2]. The 
perturbative Feynman expansion generates stranded graphs dual to simphcial complexes of all 
topologies, weighted by spin foam amplitudes [3]. Conversely, it can be shown [3] that any spin 
foam model admits a GFT formulation, which removes its dependence on the triangulation. 

Most spin foam and GFT models for quantum gravity are based on modifications of the 
Ooguri model for 4d BF theorjQ. This approach is motivated by the fact that classical 4d 
gravity can be expressed as a constrained BF theory (Plebanski formulation) [6l [7] 

S{uj,B,X)= I Tr B A F{uj) + XC{B), (1) 

Jm 

for so (4) valueci 1-form connection uj and 2-form field B, where C{B) are polynomial (so-called 
simplicity) constraints and A is some Lagrange multiplier. The variation with respect to the 



^with some exceptions, see for e.g [5]. 

^we will restrict to the Riemanian signature in this paper. 
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Lagrange multiplier constrains to be a function of a tetrad 1-form field B = *{e A e), turning 
BF to the Palatini action for gravity in the first order formalism. The Immirzi parameter 7, 
which plays a crucial role in loop gravity, can be introduced by replacing B^B + ^*B in the BF 
term of the action. Solving the constraints reproduces the Hoist action [8] , classically equivalent 
to Palatini gravity and starting point for the quantization leading to loop quantum gravity. 

The spin foam quantization stems from a discretization of the classical theory, by choosing 
a triangulation A on A^. While the most direct route to quantization would be to include a 
discrete analogue of the constraints C{B) into the definition of the measure of the discretized 
path integral [U \W\ 



the standard spin foam strategy consists of quantizing first the topological BF part of the 
discretized theory: the discretization and quantization of BF theories in any dimension are in 
fact well-understood [101 [TT] . The task is then to implement a quantum version of the constraints 
in order to recover the gravity degrees of freedom. This has shown to be a quite subtle task, partly 
because of the very simplicial setting in which the construction takes place: in fact, no standard 
canonical quantization procedure exists in such discrete setting, and of course things only become 
more difficult when the classical system to be quantized is a background independent simplicial 
gravity theorj0. Proposals for the implementation of the constraints motivated by the geometric 
quantization of simplicial structures [131 [T4] first led to the famous Barrett-Crane model [15] , 
and more recently to the EPRL model [16], which includes the Immirzi parameter and reduces 
to Barrett-Crane when 7 = 00. 

One of the main difficulties that this strategy encounters stems from the non-commutative 
nature of the geometrical variables in the BF Ooguri model [13] - beginning with the quantum 
B variables themselves represented as generators of the gauge group - which obscures the geo- 
metrical interpretation of the constraints. As an attempt to remedy this problem, Livine and 
Speziale suggested [TTj to rewrite the BF amplitudes in a basis of Perelomov group coherent 
states and to interpret the coherent state labels as classical bivectors (though with quantized 
norm) on which to impose the constraints. The realization of this idea led to the FK^ model 
[18j which, remarkably, coincides with EPRL when 7 < 1. Using the coherent sate representa- 
tion, it can be shown that both models exhibit a path integral- like formulation [19j . Though it 
provides a powerful tool to relate the models to Regge gravity in the semi-classical limit [20j . 
this formulation, which involves a quite specific and non-standard action, is however far remote 
from the original path integral ([2]). 

The goal of the series of papers [2HI22j and the present work, in the spirit of earlier works [23] 
by Bonzom and Livine, is instead to dig deeper into the relation between spin foam models and 
simplicial path integrals of the type ([2]). Our main result is to show that generic spin foam models 
based on quantum BF theory have a dual formulation as a version of ([2|) in which functionals 
of the discrete bivectors B/\ are endowed with a non-commutative structure (star product) 
deforming the usual point-wise product. Such a formulation of the path integral!, in contrast 
to standard constructions with commutative variables [23], captures the key aspect of non- 
commutativity of bivectors in spin foam models, covariant counterpart of the non-commutativity 

^see however [T^] for a recent proposal of a general canonical formalism for simplicial gravity, 
^this formulation is adapted to the quantization of classical systems with a 'curved' phase space, see [25) for 
the simple example of a quantum system on the group SO(3). 
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of the flux variables in loop gravity |26^ [27] . It is important to note that this non-commutativity 
is not an anomaly of quantization, as it can be traced back to the classical theory |26[ l28 [[29ll30j . 

This result is important for the study of spin foam models, in many respects. It gives a 
direct correspondence between purely algebraic amplitudes describing the quantum geometry 
on one hand, and a measure on the variables of the classical theory on the other. It allows for a 
direct comparison between the spin foam quantization and a proper path integral quantization 
of Holst-Plebanski gravity [311 [32] . It also opens the way for a precise study of the consequences 
of non-commutativity of the geometry inherent to these models. 

The duality between spin foam models and the path integrals ([2]) is realized at the level of 
the generating group held theories. We will thus work in this very general setting, though the 
construction could also be carried out directly at the level of the amplitudes. The mechanism 
is the following. As we will review in Section [21 in addition to its usual formulations in terms 
of gauge invariant group fields figi, ■ ■ ■ (74) or of its Peter- Weyl tensor components, the Ooguri 
GFT model for SO (4) BF theory has a dual formulation in terms of fields on four copies of the 
Lie algebra so (4), obtained by a Fourier transform: 

${xu ■■■x,):=j [d<7.]' (^(51, • • • 54) e^^^^5^ • • • e^^^^54 (3) 

which endows the space of fields with a non-commutative ^-product (dual to group convolution) . 
It can be shown that the gauge invariance of field translates into a closure condition xi + - ■ ■+X4^ = 
for the so (4) variables, which have a direct interpretation as the discrete B variables labeling 
the faces of a tetrahedron. In fact, in this representation, the GFT Feynman amplitudes are 
simplicial BF path integrals |21] . In such a representation, where all the geometrical variables 
are explicit, constrained models for gravity take a very suggestive form in terms of constrained 
fields {S -k (p)[xj) for some functions S{xj) constraining the bivectors. By construction, the 
Feynman amplitudes are simplicial path integrals for constrained BF theories. 

In principle, every GFT and spin foam model for gravity based on BF theory, thus including 
EPRL/FK, can be formulated this way, with more or less natural forms for the constraint 
functions S. In this paper, however, we rather follow a constructive approach: in Section [3l we 
define S in the most natural way in this framework: namely in terms (non-commutative) Dirac 
distributions S{xj) = 5{C{xj)) where C{xj) are the discrete simplicity constraints. The use of 
Dirac distributions effectively amounts to constraining the measure on the bivector variables. As 
we will see, the fact that the constraints are imposed on the group field, hence in all tetrahedra, 
will automatically lead to additional constraints on the connection in the path integral form of 
the amplitudes [32] . 

We will work with the linear form of the discrete simplicity constraints j33[ \TE[ [18] [34] , with 
Immirzi parameter, and a minimal extension of the group field formalism to include the normals 
to tetrahedra as an additional variable of the field: this allows us to implement the constraints in 
a gauge covariant way. In this formalism, polynomial boundary observables are labelled by so- 
called projected spin- networks [35]. We thus obtain a constrained GFT formulated as a theory 
of dynamical (non-commutative) geometric tetrahedra, which interact in the simplest possible 
way, as dictated by the star product. Its Feynman amplitudes define a spin foam model for 
gravity with Immirzi parameter 7, which gives a variant the Barrett-Crane model when 7 = 00 
but differs from the existing models for generic values of 7. In particular it does not require any 
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rationality condition for 7. This model is formulated both as a path integral ([2]) and in terms 
of Wigner 15j-symbols, in Equ. (jSip and (j52p below: this is the main result of the paper. 

The framework will also allow us to take a new (covariant) look at the peculiar features of the 
path integral amplitudes induced by non-commutativity. In particular we briefly discuss Section 
m how the so-called 'ultralocality problem' anticipated for the Barrett-Crane model manifests 
itself in our framework and argue that it may disappear in a suitable semi-classical limit involving 
a commutative limit. 

In Section [5] we conclude and sketch some directions for future work. 

2 GFT models for BF theory 

In this section, we start by recalling the standard Ooguri GFT for BF theory and its non- 
commutative bivector representation. We then present an extension of the GFT formalism, 
where the usual field variables, associated to the four triangles of a tetrahedron, are supplemented 
by an vector playing the role of the normal to the tetrahedron. As we will see, it will allow 
us to implement the linear simplicity constraints (j27p in a covariant way. 

Our notations and conventions are as follows. We identify functions on SO (4) with functions 
on SU(2)"xSU(2)+/Z2 and denote by g = {g~,g^) the SU(2) decomposition of the field variables. 
We also use the decomposition of so(4) in anti-sef dual and self dual sectors so(4) =su(2)+©su(2)" 
and denote by x = (x~,x+) the corresponding decomposition of its elements. From Sec. 12.21 on. 
based on the S0(3) Fourier transfornu [36], we further assume an invariance of group functions 
under g — )• —g, so that they are effectively functions on S0(3) xS0(3). 

2.1 Connection and spin formulations 

In the standard connection formulation, the Ooguri GFT model [TT] for BF theory is described 
in terms of a field (p{gi, ■ ■ ■ 54) on four copies of the gauge group, satisfying the gauge invariance 
condition: 

V/ieS0(4), ip{gi,---g4) = ^{hgi,---hgi) (4) 
The dynamics is governed by the action: 

5'=- y"[d5'j]'^¥'l234 + ^ y [dfifj] ^"¥'1234 ¥'4567 ¥^7389 9^96210 V^lO 851 (5) 

where 971234 is a shorthand notation for ip{gi, ■ ■ ■ (74), dg is the normalized Haar measure and A 
is a coupling constant. The perturbative expansion in A generates 4-stranded graphs dual to 4d 
simplicial complexes (see Fig[T|): if one associates the field variables to the four triangles of a 
tetrahedron, the quintic interaction sticks five tetrahedra together a common triangle to form a 
4-simplex; the kinetic term dictates the gluing rules for 4-simplices along tetrahedra. 

By using the harmonic analysis on SO (4), the gauge invariant field is expanded into four 
SO (4) irreducible representations, labelled by pairs of SU(2) spins J = and 4-valent 

intertwiners t = {i~ ,1^) labelled by a pair of intermediate SU(2) spins. The interaction vertex 

^An extension of the group transform to the whole SU(2) has been developed in [37]. We do not use it in this 
work, because on the one hand we do not expect the results to be very much modified by such extension, and on 
the other hand the general case would entail a more involved notation. Note also that a different SU(2) transform 
has been proposed and studied in |38) . 
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Figure 1: GFT propagator and vertex 



is expressed in terms of the SU(2) Wigner 15j-symbols [11] . In this formulation, the Feynman 
ampUtudes takes the form of a state sum model on the simplicial complex dual A to the graph: 

E n^.7^.-n{i5j-};{i5j}: (6) 

{JtM * 

The sum is over the SO (4) representations Jt and intertwiners ir labeling the triangles and the 
tetrahedra; dj = 2j + 1 is the dimension of the SU(2) representation j. The symbol for the vertex 
amplitude (associated to each 4-simplex), is the product of two SU(2) Wigner 15j-symbols. 

2.2 Non-commutative Fourier transform and bivector formulation 

The simplicial geometry encoded in the model ([5]) is best understood in a dual formulation, 
coined 'metric representation' in [21], obtained by a group Fourier transform of the field. The 
relevant Fourier transform here is the obvious extension of the non-commutative SO (3) Fourier 
transform [MlEZl ?] to the group [S0(3) x S0(3)]^: 

${xi, • • • X4) [dgif ifigi,- ■ ■ 54) e^^-isi . . . e^Trx«4 (7) 

The variables Xi belong to the Lie algebra so (4) = su(2) ® su(2). The kernel of the Fourier 
transform is a product of 'plane waves' Eg{x) = e*'^''^^, where the trace Tr is defined in terms of 
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the usual trace of 2 x 2 matrice^ as TTxg = Yl± eg±ti'[x^(7^] with e„± =sign(tr(7^). Thus Eg{x) 



is itself a product of two S0(3) plane waves eg±{x^) :=e*'^!'*^'"^ ^ . The plane waves satisfy the 
properties: 



I d^'xEgix) = 6{g), Eg-i{x) = Eg{-x) (8) 

Here d^x is the Lebesgue measure on so(4) and 5{g) :='5so(3)(fl'^)<^SO{3) io ) acts as the delta 
distribution on group fields. We deduce from these the following expression of the GFT action 
([5]) in terms of the dual field (p: 

•S" = ^ y[d'^3;i]^^1234 * ^-1-2-3-4 

+ ^ j [d^Xj]"^" ^1234 * <?-4567 * ^-7-389 * ^-9-6-2 10 * "^-lO-S-S-l (9) 

Notations are as follows. ^±i±2±3±4 is a shorthand notation for ^(ztxi,-- - ,ibx4). The 
product is defined on S0(3) plane waves as eg -k eg'{x) = eggi , extended to Eg -k Eg'{x) = Egg'{x) 
and by linearity to the algebra functions. In the expression above, it is understood that the 
T*r-product pairs repeated indices: for example, the first product of the interaction term is a 
product of functions of the variable X4, (^1234 *X4, ^-4567- To recover ([5]) from ([9|), one expands 
the dual fields in group modes; the integration over the variables Xj produces delta functions 
6{g~i^g'i)=J dxiEg-ig, {xi) which identify the group elements associated to the same index. 

Gauge invariance (j44p translates into the invariance of the dual field under ^-multiplication 
by a product of four plane waves Eh{xi) ■ ■ ■ Eh{xi) = Eh{xi + .. + X4) labelled by the same h: 

V/iGS0(4), ^ = Eh---Eh*^ (10) 

Integrating over h on both sides of this equality gives: 

= 5{xi+ .. +X4^) -k (f, 6{x) := J dhEh{x) (11) 

where 5 plays the role of non-commutative delta function on algebra functions 5-k(f)[x) = 0(O)5(x) . 
In words, gauge invariance corresponds to a constraint on the dual fields imposing the closure 
xi + .. -|- X4 = of its variable^. It is interesting to note that, in the non-commutative setting, 
the closure constraint is implemented by a projector, since 6 k 5 = 5. Geometrically, (p represents 
a tetrahedron whose four faces are labeled by a bivector xf'^. The glueing rules for tetrahedra 
dictated by the action ([9]) corresponds to the identification of the face bivectors, modulo a sign 
encoding a flip of the face orientation. 

More precisely, propagator and vertex in this representation are given by 

4 „ 10 

P(x,x') = J{5-x^{x[), V{x,x') = / [dh,f .J * ^h,/.-,0(^f ) (12) 

i=l i=l 



^Let Tj be i times the Pauli matrices, then txTiTj — —Sij . Given and SU(2) element u = e*"'^J parametrized by 
the angle 6 G [0, tt] and the unit R'^-vector n and a = a-'rj in the algebra su(2), we thus have tr[ait] = — sin ■ a. 
Also := sign(tru) = sign(cos 0) . 

^In terms of the canonical theory, this is indeed just the standard Gauss law corresponding to S0(4) gauge 
invariance in flux variables [27] 
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where i labels the oriented strands (triangles) and I the half lines (tetrahedra) of the graphs 
in Fig [1] and dx{y) '■= S{x — y), with 5 defined as in (j22p . In terms of simplicial geometry, 
the vertex function encodes the identification, for each of the ten triangles i of a 4-simplex, of 
the two bivectors xf, xf associated to it, corresponding to the two tetrahedra £, i' sharing the 
triangle, up to parallel transport hihj,^ from one tetrahedron to another [21] • The sign difference 
reflects the fact, in an oriented 4-simplex, a triangle inherit opposite orientations from the two 
tetrahedra sharing it. The integration over hi implements the gauge invariance (jlOp . We have 
chosen to gauge-average the vertex here, but since gauge averaging is a projection C^ = C, one 
could instead gauge-average the propagator or both vertex and propagator without affecting the 
amplitudes. 

The Feynman amplitudes are obtained by taking the ^-product of propagator and vertex 
functions, following the strands of the graph [21]. The structure of the ^-product gives a clear 
geometrical meaning to the algebraic expressions. In particular, commutation with a plane wave 
signifies a change of frame: Eh-k Lp = ^p^ -k Eh, with ip^ (x) = ip{h~^ xh) . 

For a given closed graph dual to a simplicial complex, this results in integrals over group 
elements h-j-a labelled by adjacent pairs {tetrahedron, 4-simplex}, interpreted as parallel trans- 
port from the center of the 4-simplex a to the center of its boundary tetrahedron r, and over 
so(4) variables x^jX^, interpreted as the same bivector of t seen in different frames associated 
to the tetrahedra and 4-simplices sharing t. In what follows we set /ig-r := ^ro- denote by 
h-rr' =hrahaT' the holouomy between two neighboring tetrahedra through an adjacent 4-simplex. 
The integrand factorizes into contributions of each loop of strands of the graph, dual to a tri- 
angle, taking the form of a ^-product of delta functions identifying all variables xJjX^ labelled 
but the same t, up to parallel transport between the corresponding frames. 

After integration over all variables but one per triangle xt:=x^° associated to a 'reference' 
tetrahedron ro(t), the amplitude reads: 

<Ta> t 

where Ht = /iron " " " hrj^^ro is the holonomy along the loop of A^t + l tetrahedra sharing t calculated 
for a choice of orientation and reference tetrahedror0. The integrand is the exponential of the 
discrete BF action, resulting from from the star product of Nt plane waves for each t: 

Eh^^^^^k...kEh^^^Jxt) = e^^^^^^^ (14) 

The GFT amplitudes in the bivector representation thus take the form of simplicial path integrals 
for BF theory, where field variables a;€so(4) and group elements /i€S0(4) arising from gauge 
invariance play the respective roles of discrete B field and discrete connection. 

The bivector formulation of GFT suggests clear routes for defining geometrical models, by 
means of constraints operators implementing the simplicity constraint on the field variables. 
The linear constraints (|27p . however, involve an another geometrical variable: the normal to the 
tetrahedron. In the next section, we review an extension of the usual GFT formalism introduced 
in [22], which include the normals as an additional field variable. Although the extended GFT 
generate the same BF amplitudes for closed graphs, it will allow us to impose the constraints 
covariantly, that is on fields that are gauge-invariant under a simultaneous SO (4) rotation of 
both bivectors and normal vector. 

^The amplitude does not depend on these choices. 
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2.3 Introducing normals: extended GFT formalism 

In the extended GFT formalism |39tl22j. the basic group field '^kidi-,' ' ' Sa) is supplemented with 
a fifth variable k G SU(2) 5^, viewed as a unit vector in M^. In geometrical models, k will 
be interpreted as the normal to a tetrahedron. Gauge invariance (j44p is replaced by a gauge 
covariance with respect to the normal k: 

Mh, ^Pk{gi,- ■ ■ Oi) = V'h^kihgi,- ■ ■ hgi) (15) 

where h> k := h+k{h~y^ is the normal rotated by h. Clearly, the field obtained by integrating 
over the normals obeys the gauge invariance (j44p . 
The dynamics is governed by the action: 

S[V>] = \ j [dgifdk^pl^2M + [d5i]^°[d^i]^V'fcil234<^fe24567¥'fe373899'fe496210V'fc5l0851 (16) 

where V'fci234 is a shorthand notation for 97fc(5'i, • • • (74), dg and dk are the Haar measures on 
S0(4) and SU(2). Hence, whereas the interaction polynomial does not couple the normals, 
the kinetic term, which encodes the glueing rule of 4-simplices along a tetrahedron, identifies 
both group elements and normals. It is already clear from the structure of this action that the 
amplitudes of closed diagrams will not depend on the normals; the extended formulation only 
modifies the structure of boundary states. 

Note that gauge covariance (jl5p induces an invariance under the stabilizer group 80(3)^ = 
{/i€ S0(4), /i+A;(/i~)~^ = /c} of the normal /c, affecting only the four group arguments of the field. 
Upon Peter- Weyl decomposition, gauge invariant fields are expanded into four irreducible SO (4) 
representations (given by pairs of SU(2) spins Jj = (j7,ij^), i = 1 ■ ■ ■ 4), each of which can be 
further decomposed into 80(3)^ representations. A set of basis functions is given by: 

\t=l / 

where repeated lower indices are summed over. (5±) are the 8U(2) Wi gner matrices, (ij) 
form a basis of four-valent 80(3) intertwiners, labelled by an intermediate spin j. The k- 

dependent coefficients, defined in terms of the 80(3) Clebsch-Gordan coefficients CtnLp^ as: 

m 

define a tensor that intertwines the action of 80 (3) ^ in the representation and the action 

of 80(3) in the representation ki. Namely, given Uk = {k~^uk,u)GS0{3)k, we have: 

C'l^'^^J ik)D'\ _(u,)D^V +(u+) = C^I^v''(A:)L»„%.(n). (19) 

()17p corresponds to the vertex structure of the so-called projected spin networks of the covariant 
approach to loop quantum gravity [35], which thus define a basis for polynomial gauge invariant 
observables (and thus boundary states) in the extended GFT formalism. 

Just as in the standard formulation, the bivector representation of the GFT is obtained by 
Fourier transform of the field 

ipk{xi,..,x^):= J [d<7]Vfc(5i,-, 54) 6^1^x454 (20) 
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Gauge invariance translates into: 



V/i, = Eh - ■ ■ Eh* ^h-i>fc (21) 

where h~^>k = {h^y^kh~ , thus implemented by the gauge invariance projector acting on extended 
fields as 

{C>lp)k = j dhEh---Eh-k (ph-i^k (22) 

Note that upon integration over the normal, gauge invariance gives the closure of the four 
bivector variables: in fact if = / dkipki then ■0 = 6{xi + .. + X4) -k'll). The action is the obvious 
extension of ([9]). It is interesting to write the interaction in terms of ■0, to emphasize the fact it 
implements the closure constraints: 

^ = \ / ^fcl234 *^fc-l-2-3-4 
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+ ^ / [d''Xi]^°'i/'1234 *V'-4567 * V'- 7-389 *V'-9-6-2 10* fA-lO-S-S-l (23) 



The propagator of the extended GFT is supplemented with an additional strand which 
identifies the normals k up to parallel transport arising from gauge invariance. Just as in the 
non-extended case, gauge invariance can be implemented in the vertex only, or in the propagator, 
or in both vertex and propagator. Choosing the first case, the propagator then reads, in the 
bivector formulation: 

4 

P{x,k-x\k') = \{5-^M)Kk'k-^) (24) 

i=l 

The additional contribution reduces the number of SU(2) variables to one per link, hence, in 
terms of the dual simplicial complex, to one k^ per tetrahedron. The integrals over the normals 
on the internal links (bulk tetrahedra) drop from the amplitudes; just like in the standard 
formulation, the extended GFT generate simplicial BF path integrals as Feynman amplitudes. 
However, due to the extension of the GFT field to include the normal vectors to tetrahedra, 
the boundary states appearing in the amplitudes for GFT n-point functions are differerent from 
those of the standard Ooguri model, recovered only after averaging out the normal variables 
independently at each tetrahedron in the boundary. 



3 Geometrical models 

We have seen in the previous section that, in the bivector representation of the GFT for BF 
theory, the field ip = J dkipk represents a tetrahedron characterized by four bivectors Xj, j = 1...4 
playing the role of discrete B field; gauge invariance implements the closure condition V • x,- = 0. 
We now propose a natural modification of the GFT (I16p in terms of a constraint operator acting 
on the field by implementing the simplicity constraint of its bivector variables, for any positive 
value of the Immirzi parameter 7, allowing to reconstruct a discrete tetrad for the tetrahedron. 
By construction, the Feynman expansion will generate simplicial path integrals for a constrained 
BF theory of Holst-Plebanski type. 

We start by recalling the discrete form of the (linear) simplicity constraints for a classical 
bivector geometry. 
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3.1 Discrete simplicity constraints 

In the absence of Immirzi parameter, the simphcity constraints state that the hodge dual *xy 
are the area bivectors of a geometric (metric) tetrahedron: this is the discrete equivalent of 
B = *e A e. Following \32\ I18j . these constraints are implemented by requiring that the four 



*xy lie in the same hypersurface normal to a given unit vector in R^, namely *xykj = 
for all j = l..A. Using the selfdual/anti-selfdual decomposition of the algebra so(4), this can be 
expressed as: 

Vj G {1...4}, 3k G SU(2), kx'jk'^ + x+ = (25) 

The variable k € SU(2) ~ is then the SU(2) representation of the unit vector k^ normal to 
the tetrahedroill]- If ten bivectors labeling the faces of a 4-simplex satisfy simplicity and closure 
constraint for each tetrahedron, then they define a geometric 4-simplex (for non-degenerate 
configurations). Furthermore, if the (constrained) bivectors associated to a given tetrahedron 
are also correctly identified across the two 4-simplices sharing it, then the reconstruction of a 
discrete tetrad can be carried out for the whole simplicial complex, again modulo degenerate 
configurations. 

The inclusion of the Immirzi parameter can be performed easily also at the discrete level 
|18j . We have mentioned that the Immirzi parameter is introduced in the continuum action ([T]) 
by a change of variables B ^ B = B + ^*B in the BF term. The action in the new variables 
looks again like a constrained BF theory, but where the constraints C{B{B)) are now imposed 
on the following linear combination field B: 

B{B) = —^i*B-^B) (26) 

In the discrete setting, the simplicity condition with Immirzi parameter is thus obtained from 
(j25p by replacing by replacing xj by the linear combinations *xy — ^xY . In terms of the 



selfdual/anti-selfdual decomposition of the bivectors, it reads: 

Vj e {1...4}, 3k € SU(2), kx-jk-^ + /3x+ = (27) 
where the parameter /3 is related to the Immirzi parameter as: 

/3=^ (28) 
7 + 1 

Note that the relation (I27p is invariant under simultaneous sign flip 7 — >■ —7 of the Immirzi 
parameter and exchange x'j o x~- of the self dual and anti-self dual part of the bivectors. 



It is also well-defined for 7 = 1, although then the change of variables (|26p is singular and the 
geometrical interpretation is lost. In the following, we restrict to 7 G [0, 00], so that the parameter 
/3 takes its value in [—1,1]. 



3.2 Constraint operator and non-commutative tetrahedra 

Back to GFT, we now need to encode the simplicity condition (I27p of the bivector variables Xj 
as a constraint on the field (pk, the idea being of course to identify the normal to the tetrahedron 
to the SU(2) variable k of the field. The natural way to do so why taking into account the 

^Let k-—{k~,k^) be the SO(4) rotation mapping the vector A^^ = (1, 0, 0, 0) to k' , then k^k'^k'^. 
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non-commutativity of the fields is to use non-commutative delta functions, defined by their 
plane wave expansion (j22p . These delta functions act as distributions for the star product: 
6 -k (p(a) = (j){0)6{a), so using these to constrain the field will effectively amount to constrain the 
measure on the bivectors. 

We thus introduce the following function of x = {x~,x^) G 5o(4): 

S^(x) := 6_k^-k.i{Px+) = [ dne^*'^['=''"'=^"]e^^t''["^''] (29) 

^SU(2) 

where S-a{b) := 5{a + 6) and 5 is the su(2) non-commutative delta function. Our geometrical 
GFT models will be defined by constraining the field (pk{xj) in the action by means of an 
operator acting on it by ★-multiplication by the product S^{xi)...S^{x4) of four simplicity 
functions: 

4 

{S/^ >ip)k{xi, - ■ ■ ,X4^) = Y[S^{Xj)-k(pk{xi,...X4) (30) 

i=i 

We give below the explicit expression of the star product (j30p in terms of group and tensor 
Fourier components. But first, let us show that the action of this operator is well-defined. We 
will also see that it commutes with the SO (4) gauge transformations ()2ip . 

To be able to take the ★-product of the with the field, we need the function (j29|) to be 
in the image of the group Fourier transform. To see why this is indeed the case, notice that, 
because |/3| < 1, there exists u^€SU(2) such that /3tr[ati] =tr[an^] for all a € su(2). Indeed, if 
u = e^'^^'^j is parametrized by the angle 6 € [0, vr] and the unit M^-vector n, where the Tj are i 
times the Pauli matrices, we define up = e^^'^^'^\ where the parameters Op and np are 



sin6'^ = |/3| sin^, sign(cos ^/j) = sign(cos ^); = sign(/3)n (31) 
The simplicity function (j29p can thus be written as a superposition of plane waves Eg{x) = e^^'^^^'. 



(/3x+)= / AuE.{x) (32) 

iSU(2) 

where we introduced = {k~^uk, up) G SU(2) x SU(2). Therefore it belongs to the image of the 
Fourier transform, and its star-product with the field is well-defined. 

The operator S is not a projector for generic values of /3, unless /3 = 0, 1 (which corresponds 
to 7 = l,oo). Indeed, because of the nonlinearity of scaling by /3 in the definition (j3ip . we have 
that {uv)p^upvp, and thus * S^=/S^. Remarkably, however, the action of S is well-defined 
on gauge invariant fields, as it commutes with the gauge transformations 

S^>[Eh---Eh* = Eh---Eh*{S^> (33) 

thanks to the commutation relations between plane waves and simplicity functions: 

Eh*S§ = S^,^,*E, (34) 

Geometrically, these relations express the fact that rotating a bi-vector which is simple with 
respect to a normal k gives a bi-vector which is simple with respect to the rotated normal 
h> k := h^k{h~)~^ . This is the advantage of the extended GFT formalism, where the normals 
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are explicit variables of the field: the linear simplicity constraints on the bivectors can be imposed 
a covariant way. This is not the case in the standard formulation of the Barrett-Crane model, 
nor on the EPRL-FK model, where simplicity and gauge invariance are implemented by means 
of two non-commuting projectors. 

Let us now examine the dual action of S on the original group fields ^kiSj)- By using the 
plane wave expansion (j32p of the simplicity functions and the definition of the star product, we 
obtain: 

{Sf^ >ip)kigi,---,9i) = [duj]''(^fc(ui«5'i,---,uL5f4) (35) 

where u*^^ = {k-^ujk, Ujfs) € SU(2) X SU(2) and Uji3 is defined as in (I3ip . For the particular value 

/3 = 1 reached in the limit 7—7-00, reduces to the projector onto fields on four copies of the 
homogeneous space SO(4)/SO(3)fc. Using the invariance (jl5p . one can gauge fix the normal to 
the value k = 1 (time gauge). On such gauge fixed fields, the simplicity operator coincides with 
the projector defining the standard GFT formulation of the Barrett-Crane model [l5l[22]. The 
diff'erence here, however, is that the gauge fixed extended fields are obviously not gauge invariant 
under the full S0(4) - but only under the diagonal SU(2) subgroup. 

Upon Peter- Weyl decomposition of the constrained field {S^ >y^)kj a set of basis functions is 
given by the action of the S on the functions (fT7|l : 



where repeated lower indices are summed over. This expression is obtained from (I17p by replacing 
the k-dependent coefficients C'' ^ ^ ik) = C^ ^ _ (k) by new ones given by: 

F^l^\^\k)= I duD^\ .{k-^uk)D^\ +{ub)C^~^^^^: {k) (37) 

with Ufi given as in pip . Just as in (|17p . these coefficients intertwine the action of stabilizer 
subgroup S0(3)jt in the representation H and the action of S0(3) in the representation ki. 
Namely, given = (A;"^mA;, u) G 80(3)^, we have: 

F^l^'^J {k)D'\ Ml)D'\ +{vil) = F'['\\k)D%Xu). (38) 

Here, they also contain all the information about the simplicity constraints and the specific form 
of the operator that implements them. In particular the integral of the two Wigner matrices 
encodes a relation between the spins which depends on the Immirzi parameter; for 

example j~ when /3 G {-1, 1}, namely when 7 € {0, 00}. However for generic values of /3, it 
does not enforce the spin relations j~ = characteristic of the EPRL-FK models, an analogue 
of which we expect to recover only in the asymptotic regime. In particular, they do not impose 
any rationality condition on the Immirzi parameter 7. A detailed study of the properties of 
these coefficients, in particular their asymptotic behaviour for large spins, is left for future work. 

To sum up this section, we have defined an operator 5^ acting on the field by imposing 
the linear simplicity condition (j27p . for any positive value of the Immirzi parameter, on its four 
bivector variables. Note that, because of gauge invariance, the closure constraint holds after 
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integration over the normal: let ■0 := J dkS^ t> ifk, then ^ = 5{xi + .. + X4) * where 6 is the 
non-commutative delta function 5 -k (p{x) = 4>{0)6{x) defined in (j22p . 

In relation to a canonical theory, on cane link the definition of the constrained GFT field 
to a quantization of a tetrahedron characterized by its four constrained bivectors: in particular, 
one can check (using the Fourier duality with group fields) that the generators Jj'^ of the gauge 
group act by multiplication by the coordinate functions xY {xj) = xY . The quantization 
procedure consists of first quantizing classical configurations {xj} G 5o(4), k G SU(2) ~ of 
bivectors and normal and then to impose geometricity (simplicity) constraints at the quantum 
level. The Hilbert space is the tensor product of L^(SU(2)) with 

4 

i=l 

where the spaces, which also appear as state spaces in the flux representation of loop quantum 
gravity [27], are spaces of functions on so(4) endowed with the scalar product J d^x{f-kg){x) 
where f{x) = f{—x). The algebra structure encoded in the star product, which deforms the usual 
point-wise product, stems directly by Fourier transform from the algebra of group functions. It 
makes explicit the non-commutativity of the geometry inherent to spin foams and group field 
theories. This procedure is manifestly dual to geometric quantization: the advantage here is 
that the classical variables characterizing the geometry remain explicit, as arguments of the 
fields. Geometricity conditions are then implemented by using two commuting operators: the 
simplicity operator and the gauge projector C defined by ([22]) . leading to the unambiguous 
definition of a 'geometricity operator' G = S^C = CS^ . All constraints are imposed by means of 
non-commutative delta-functions, acting as Dirac distributions for the star product, so that it 
effectively amounts to constrain the measures d^Xj on the classical field variables. 

The algebra structure and scalar product allow to build up more involved polyhedra obtained 
by glueing tetrahedra along a common triangle. Thus, glueing five geometric tetrahedra ^ := 
J dkG i^k along common triangles as 

$1234 * $-4567 * $-7-389 * $-9-6-210 * $-10-8-5-1 , (39) 

where \I'/;-|-i-i-2it3±4 is a shorthand notation for ^'^(ibxi, • • • ,±Xi) and the star product pairs 
repeated lower indices, gives a straightforward ansatz for a 'quantum 4-simplex'. This is precisely 
the interaction polynomial of our geometrical GFT, which we define now. 

3.3 GFT for gravity with no Immirzi parameter 

We first focus on the particular case (3 = 1. It corresponds to 7-^00, namely the case of gravity 
with no Immirzi parameter. In this case, the operator reduces to the projector onto the 
homogeneous space SO(4)/SO(3)fc. The corresponding geometrical GFT has been introduced 
and studied in [22]; we recall the basics here, before treating the general case. 

The model is defined by constraining the field (pk in the action of the extended Ooguri model. 
Since is a projector, it can be inserted in the propagator, in the vertex, or in both - without 
affecting the amplitudes. In particular, these choices all lead to a unique form of edge amplitudes 
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in the spin foam model. Here we set ^I':=/ dkS^ o (pk and define tlie action: 

^ = \ j [d*^^*]^ d/c ^fel234 * '^fc-1-2-3-4 

1234 * ^-4567 * ^-7-389 * ^'-9-6-2 10 * *-10 -8-5-1 (40) 

where the star product pairs repeated indices. 

We now calculate the Feynman amplitudes of this theory. We will obtain two dual represen- 
tations of these, in terms of simplicial path integrals on one hand, and spin foam models on the 
other. 



3.3.1 Simplicial path integral representation of the amplitudes 



The Feynman amplitudes of this theory are calculated with the same propagator as in (j24p and 
the vertex: 



^ 10 

V{xi- h) = / [d^s Yl{S_,, * Sk, * Ef^^^-j * )(xf ) 
1=1 



(41) 



where i labels the oriented strands (triangles) and £ the half lines (tetrahedra) of the vertex 
graph, and 5"^ = (5_;.^.-^-i (x+) is the simplicity function for /3 = 1. The calculation is analogous to 
the unconstrained case. For a given closed graph dual to a simplicial complex A, this results in 
integrals over holonomies hra G S0(4), tetrahedra normals kr G SU(2), and bivectors xj, x" on t 
seen in different frames associated to the tetrahedra and 4-simplices sharing t. After integration 
over all variables but one per triangle xt:~'^^° 
amplitude reads [22] : 



associated to a reference tetrahedron ro(i), the 



-dK^[<lk,][d^xt] 



*j=0 Shojt>kj {xt) 



★ e 



(42) 



The notations are that of Sec l2.2l Ht = hj 



is the holonomy along the (oriented) loop 



of A't -|- 1 tetrahedra sharing t, labelled by the integer j; h^j = h-rgn ■ ■ ■ ^rj_iTj is the holonomy 
from the reference tetrahedron to the j-th tetrahedron sharing t. kj '■=kr^ is the normal of the 
j-th tetrahedron around t. The function Sh^^.^hjixt) imposes on xt the linear simplicity condition 
with respect to the rotated normal h^j > kj := hQjkj{hQj)~^ , namely the pull back of kj in the 
frame of the reference tetrahedron To(t). It amounts to imposing the linear simplicity of the 
pushed forward bivector h'^jXthoj with respect to kj. The integrand results from taking, for each 
t, the alternate star product of Nf plane waves and iVj simplicity functions: 



Sko * E, 



hr 



^ j=0'-'hojl>kj 



* e 



iTvxtHt 



(43) 



where we used the commutation relation -k Ef^ = Sh:>k * Efi to regroup all the plane waves on 
the right of the expression. 

The Feynman amplitudes of the GFT ([9]) thus take the form of (non-commutative) sim- 
plicial path integrals for a constraint BF theory of Plebanski type. The constraints are non- 
commutative delta function modifying the measures d^Xt on the bivectors, imposing the sim- 
plicity of each xt with respect to the normals of all the tetrahedra sharing t. 
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Note that by construction, the integrand of [32] is invariant under S0(4) rotations {gnOa} of 
all local frames: 



hra ^ grhraOa, K H> §+^{9^)'^ Xf ^ 9r\t)Xtgr(t) (44) 

where T{t) is the reference tetrahedron of the triangle t. This includes the gauge invariance of 
the discrete BF action. The choice Qr = {g'r, dt)'-— iK^ > 1) leads to the 'time gauge' A;,- = 1, which 
shows that the integral over the normals /cr drop out of the amplitude. Of course in the case 
of open graphs, dual to simplicial complexes with boundary, the amplitude still has an explicit 
dependence on the normals of the boundary tetrahedra. 

It is interesting to distinguish two types of constraints on the bivector xt of a given triangle. 
The j = contribution S'fcg(a;() = 5_^^^.-^-i(a;^) in (jl2]) imposes the linear simplicity of each bivector 
Xt with respect to the normal to the reference tetrahedron To(t). The remaining part, for a given 
set of bivectors, can be viewed as constraints on the holonomies modifying the measures d/iT-o- 
on the discrete connections. The effective measure 



transforms covariantly under the gauge transformations (j44p . Using the cyclic invariance of the 
★-product under integration, the amplitude can be written in terms of this measure as: 

where the star product pairs the variables xj. Note that the appearance of a gauge covariant 
measure on the discrete connection is the result of using extended GFT fields, and thus of 
requiring a covariant imposition of the simplicity constraints. The need of a generalization of 
the closure constraint to achieve this has been noted on several occasions in the literature |321l34j 
(see also [3Tj). 

3.3.2 Spin foam representation of the amplitudes 

We have derived the GF amplitudes starting from the bivector formulation of the group field 
theory, where the simplicial geometry is implicit. By construction, they take the form of simpli- 
cial path integrals. It should be clear, however, that the dual connection and spin formulations 
of the same GFT will dual expressions of the same amplitudes in terms of a lattice gauge theory 
and spin foam amplitudes. 

The spin foam representation of the amplitudes (142p can be also be computed directly by 
Plancherel decomposition of the group functions into irreducible representations and integration 
over group and Lie algebra elements. As proved in [22], it gives the Barrett-Crane amplitudes: 

xa = tbc := e n 4 n Tr\r n^io^-}^ (47) 

The sum is over SU(2) spins jt labeled by triangles; dj^ := 2jt + 1. The products are over all 
triangles t, tetrahedra r and 4-simplices a. The 4-simplex weight {lOjjo- is the Barrett-Crane 
lOj-symbol [15]. This derivation singles out a specific edge amplitude (tetrahedral weight), which 
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differs from the ones that appear in the hterature with the exception of [23]. We emphasize 
again that, in our GFT construction, although the amphtudes of closed graphs reproduce the 
Barrett-Crane amplitudes, the requirement of a covariant imposition of the simplicity constraints 
imposed the use of extended boundary states which include tetrahedra normals, hence labelled 
by projected spin networks. 

This geometrical construction sheds an interesting new light on the Barrett-Crane model. In 
particular, the formula (I42p gives a new simplicial path integral formulation of its amplitudes, 
making the simplicial geometry of the model manifest. We refer to [22] for a detailed analysis 
and discussion. 



3.4 Including the Immirzi parameter: GFT for Holst-Plebanski gravity 

We now turn to the inclusion of the Immirzi parameter in the model. For generic values of 
13 = Ip^, the operator is no longer a projector. Depending on whether it is inserted in the 
propagator, in the vertex, or in both the vertex and propagator, in a single of in multiple copies, 
will lead to a priori different spin foam amplitudes. These will however have the same vertex 
amplitude (4-simplex weight) and differ only in the edge and face amplitudes (weights associated 
to tetrahedra and triangles). Here, just as in the previous section, we constrain the field in the 
interaction of the extended Ooguri model. Setting ^'^ := J dA;5^ > ^fcj we thus consider: 

1234 * ^-4567 * ^-7-389 * ^-9-6-2 10 * ^-10-8-5-1 (48) 

where the star product pairs repeated indices. 



3.4.1 Simplicial path integral representation of the amplitudes 

The derivation of the Feynman amplitudes is analogous to the case /3 = 1. They are calculated 
with the same propagator as in (j24p and the vertex: 

„ 10 

Vf'ixh ke) = / [dhef nC'^-xf * ^£ * Eh,h-} * Si){xt) (49) 
i=i t I 

where i labels the oriented strands (triangles) and i the half lines (tetrahedra) of the vertex 
graph, and = 6_i^^-i,-i{f3x^) is the simplicity function. Using the same notations as in (^2]) . 
we obtain, for the amplitude of a closed graph: 

(50) 

where S^*^ denotes the squared function 5^ * 5^. For each triangle t, the constraints impose, 
by means of non-commutative delta functions, the linear simplicity condition of Xt with respect 
to the normals of all the tetrahedra {Tj}j=o,„N^ sharing t. The square stems from the fact 
both 4-simplices sharing the tetrahedron tj contributes to a factor S^^^^j^Xxt). The Feynman 
amplitudes of this theory thus take the form of simplicial path integrals for a constrained BF 
theory of Holst-Plebanski type with Immirzi parameter 7, with linear simplicity constraints [34] • 



n 



Nt ql3*2 

j=0 "^hojt-kj 
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Just as for the case with no Immirzi parameter, the integrand of (jSOp are invariant under 
the gauge transformations (j44p : in the case of closed graphs, the integration over the normals 
k-r € SU(2) drops from the amplitude. Making explicit the form of the simplicity functions, we 
thus obtain: 

where we wrote for the squared function / * /, and hoj = h^-{hQ^)~^ . Finally, by splitting 
the constraints into a part (the i = contributions) that is independent of the holonomies and 
a part playing the role of constrains on the holonomies, the amplitude z£ can be put under a 
form analogous to (j47p . in terms of a covariant measure on the space of discrete connections. 

3.4.2 Spin foam representation of the amplitudes 

The spin foam representation of the amplitudes can be obtained either directly from (j50p by 
inverse Fourier transform and Peter-Weyl decomposition of the group functions, or from the 
Feynman rules of the spin representation of the generating GFT. In this section we give the 
explicit form of the resulting spin foam amplitudes in terms of 15j-symbols and so-called fusion 
coefficients, which will allow a direct comparison with the existing models [16[ [TSt 119] . Their 
derivation from the GFT is straightforward: we only sketch it here. 
In the spin foam representation, the amplitudes read: 

^A= E J\d^-d.+ \{d^,^\{Al{jt,Ki,-K) (52) 

where the 4-simplex weight (vertex amplitude) is given by: 

AP{jtAtr,ir; K) = E {15i};{15i}: n d.,^^d^+Jl: + (j±, Xtr; K) (53) 

The notations are as follows, t, r, a denote the triangles, tetrahedra and 4-simplices of the simpli- 
cial complex A. The sums are over S0(3) representation j, A and four-valent S0(3) intertwiners 
I, all labelled by an integer spin, and dj = 2j + 1. This gives a pair of spins {jt,jf) for each 
triangle, a spin Afr for each couple (tr) with t C r, a a spin i-j-o- for each tetrahedron and a pair 
of spins for each couple (ra) with t C a. We set dj = 2j + 1. The variables kr G SU(2) 

are the normals to the tetrahedra: as we have seen, the dependence upon the normals for the 
bulk tetrahedra (internal links of the GFT graph) drop, hence we haven made the integrals over 
these explicit. The above amplitude may as well be evaluated in the time gauge /c,- = 1 for all 
bulk tetrahedra, though it then makes less transparent the nature the boundary states, here 
labelled by projected spin networks. 

The amplitude is defined in terms of SU(2) Wigner symbols {15j}^{j^,i^^) and so- 
called fusion coefficients jl6j /*J ^ . These coefficients define a map from the space of S0(3) 

intertwiners between the representations Xtjr, ■■, ^Ur and the space of S0(4) intertwiners between 
the representations {Jl 



[d/i..][d6xt] 
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While the form (|52p is quite general for a spin foam model defined a a constrained BF theory, 
the specificity of a model lies into the exact form of the fusion coefficients, which encode the 
way simplicity constraints are imposed. For the new model presented here, they are given by: 

= (v.L-(^^.w(ni^;S;r(^.))(vk (55) 

where repeated lower indices are summed over. The tensor F■'^^■'^^^^^ {K) are the ones defined 
in (j37p . They provide an embedding of SO (3) structures into SO (4) ones. In particular, because 
of the intertwining property ([38]) . they realize 80(3) as the stabilizer subgroup SO(3)fcCSO(4) 
of the normal kj- to the tetrahedron. They also depend on the Immirzi parameter and encode 
the simplicity constraints. 

This form (|52p of the amplitudes follow from the GFT Feynman rules in the spin represen- 
tation. In this representation, the bivectors xt are replaced by pairs of spins Jt ■= [jl^Jt) ^'^^ 
magnetic numbers labeling the strands of the graphs. A way to read these rules is then the 
following. For a given graph Q labelled by spins Jt and tetrahedron normals kr, they attach 
to each 4-stranded line (tetrahedron) a propagator P^-^i{kr) G End((g)^^i JtJ defined as an 
endomorphism of the tensor product of the representations labeling its strands (triangles). The 
amplitudes are obtained by taking the trace (ie. index contractions) of all propagators, following 
the combinatorics of the graph and by summing over all spins (and normals) as: 



^A = EnV.-Tr, 

{Jt} * 

where dj = 2j + 1. The propagator decomposes as 



PP^-^'{kr) 



(56) 



pP"'^{kr) = Ps6i4)P^'''if')Pshi4) (57) 

where -fgo{4) projector onto S0(4)-invariant tensors and P^''^^{k) is defined in terms of 

the tensors (f37|) as: 

^t-n:,.^(^) = n [T.'^^^,.(^)^!^f:s^ ] (Ps%isM (58) 

1=1 \ Xi 

where repeated lower indices are summed over. -Pgo(3) is the projector onto S0(3) invariant 
tensors in Note that its insertion is actually redundant in the definition of the propagator: 

indeed, because of the property ([55]) . an SO (3) rotation in ®iXi is intertwined by F with a 
S0(3)fc rotation in ®iJi, which can be reabsorbed into P^Q^^^y This form however allows us to 
split the trace in (j56p into a product of 4-simplex weights as in (j52p . This is done by expanding 
the projectors -fgo(4) ^"^^ -^so(3) valent intertwiners {I'^i'^) and i. 

To close this section, we emphasize again that, as in the extended BF case, the boundary 
states of the model are by construction (constrained) projected spin networks [35]. Even in 
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the presence of the Immirzi parameter, we see therefore that the boundary states of the amph- 
tudes (or the GFT polynomial observables) are different from the states of standard LQG. The 
projected spin network structure is actually present also in the boundary states of all the new 
models |16^ [18] . even if in a less explicit way, and even if their apparent coincidence with LQG 
states in representation space (due to the specific form of the simplicity constraints imposed 
there) is more emphasized. 

3.5 Limiting cases 

The model presented in the previous sections corresponds to a candidate quantization of a 
simplicial version of the Plebanski-Holst formulation of 4d gravity, for generic values of the 
Immirzi parameter 7. We have already seen that the case f3 = l, which corresponds to 7 = 00, 
gives a variant of the Barrett-Crane model with a specific edge amplitude, where the boundary 
states are extended to include tetrahedra normals, hence are labelled by projected spin networks. 
Thus, the formula (j42p not only gives a new simplicial path integral formulation of the BC model, 
but it also provides a natural deformation of that model which includes the Immirzi parameter. 
We now discuss briefly other limiting cases of this model: the 'self-dual' case 7=1 and the 
'topological' case 7 = 0. 

For 7 = 1, as mentioned the change of variables [26l becomes singular - so the contact with the 
classical Hoist theory is lost. Despite the lack of a clear geometric interpretation, the constraint 
operator and the resulting model are well-defined for (3 = 0; it acts on the field by (pk by 
projection of its bivector variables onto the selfdual part of so (4). The constrained model reduces 
to the Ooguri model for topological SU(2) BF theory. 

The case 7 = corresponds to the so-called topological sector of Hoist gravity. This denomi- 
nation comes from the fact that the term of the classical Hoist action that seemingly dominates 
in this limit is the one that vanishes on shell, due to the requirement of torsion free-ness of the 
connection. As a consequence, one would expect that the resulting spin foam model/path inte- 
gral would define a trivial dynamics for any boundary state. It is not totally obvious, however, 
that the above reasoning goes through in the quantum theory as well. It could also be argued 
|42j that the resulting quantum theory would rather correspond to a quantization of 2nd order, 
metric gravity with no torsion. The rough argument is that in a path integral for the Hoist 
action, the limit 7 — )• would force, analogously to a semi-classical limit, the same path integral 
to be dominated by solutions of the equations of motion coming from the "topological term" only, 
that is exactly the torsion free-ness condition. While these arguments are obviously not con- 
clusive, they suggest not to dismiss the resulting model as un-interesting. In our context, this 
corresponds to (3 = —l. As discussed also in (22], the constraint operator for /3 = — 1: a) projects 
onto simple S0(4) representation J = (j, j), b) does not impose any restriction on the expansion 
of (j, j) into SU(2) irreducible representations A; = 0, • • • 2j, and acts on each component {J,k) 
by multiplication by the phase (—1)^-'"'"'^. In computing the amplitudes, the phase factors cancel 
each other. The resulting model, distinct from the EPR amplitudes [16], is obtained from the 
S0(4) Ooguri model ([6|) by restriction the representations to simple ones Jt = {jt, jt)- It would 
be interesting to study what the geometric interpretation of such amplitudes may be. 

The general GFT model for arbitrary f3 thus encompasses and generalizes several distinct 
models, and interpolates between them. It is tempting to speculate (see also [1I1112]), that the 
model possesses a non-trivial renormalization group flow in parameter space (Gat, /?), where Gtv 
is the Newton's constant, which is hidden in our formulation as we use dimensionless quantities 
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throughout. The natural candidates for fixed points would then be these special values for 
the (3 parameter: (—1,0,1), namely 7 = 0, 1 or 00. In particular, while the case /3 = —1 is 
distinguished only for being in some sense 'extremal', and for its peculiar classical analogue, 
the other two values can be seen as special already at the level of the very definition of the 
corresponding model. In fact, as we have seen, in these two cases, and only then, the simplicity 
operator defines a projector, and the quantum amplitudes are insensitive to the specific choice 
of insertion of this operator in the GFT action. 

4 Simplex correlations and ultralocality 

The main advantage of the framework developed in this paper is that the variables encoding the 
(fuzzy) bivector simplicial geometry of GFT and spin foam models are explicit. In particular, 
with respect to other formalisms, it gives a more direct access to the way simplices are correlated 
in the model, namely, how the model relates the geometrical data of common subsimplices in 
the glueing of neighboring simplices. 

It has in fact often been argued that, for example, the Barrett-Crane model, suffers from 
a default of correlations between neighboring simplices. This 'ultralocality' feature has been 
one of the reasons to discard this model in favour of the new models. Since the Barrett-Crane 
amplitudes show up in our model for the value 7 = 00 of the Immirzi parameter, this feature 
can in principle be clarified and dwelled further. An indepth study of the ultralocality issue is 
beyond the scope of the present paper. In this section we however discuss how it manifests itself 
in our framework. 

The interaction and kinetic polynomials of our GFT model are written as a simple star 
product of copies of the constrained field in which the bivector variables associated to the 
common triangles are strictly identified (modulo an orientation flip). At the level of the Feynman 
amplitudes however, after expansion of the gauge invariance operator, the Feynman rules dictates 
the relation between bivectors {xj",^^} expressed in different frames related by holonomies /iro-- 
For example in the BF model, two bivectors xj ,xj on the same triangle t but seen from different 
tetrahedra are related by a non-commutative delta function: 

{S-,r^Eh^^,){xl') (59) 

where /i^,-/ = hT-ijhaT' parallel transports the frame of one tetrahedron to the frame of the other. 
The structure of the star product gives a clear geometrical meaning to the algebraic expressions. 
In particular, the algebraic operation corresponding to the parallel transport of bivectors is the 
commutation with plane waves: Eh * f = f'^ * Eh, where f^{x) = f{h~^xh). 

In the amplitudes of the constrained theory, the plane waves are supplemented with simplicity 
functions imposing the linear simplicity condition of the bivectors in each frame: 

Eh^^> Sl{xt)^Eh^^,*Sl^{xt) (60) 

so that the commutation with the plane wave encodes the parallel transport of simple bivectors. 
In the case /? = 1 (ie 7 = 00) corresponding to the Barrett-Crane amplitudes, the simplicity 
functions Sk-=Sl. satisfy Sk * Sk = Sk- The definition (|29l) of these functions and the structure 
of the star product then lead to the identity: 

Sk^ * Eh^^, * Sk^, = Sk^ -k E^rf^^^^^T' * Sk^, (61) 



20 



for all uj = {k~^ulkr,ul) and u[ = {k'jiuj kr',uj ) in the stabilizer subgroups S0(3)a:^ and 
S0(3)fe^, of the normals. This is because u[,uj' can be reabsorbed into the group elements 
labeling the plane wave expansion (j29p of S^^ and S^^,- This identity can be understood as 
a relaxation of the parallel transport condition, or a weakening of the bivector correlations: 
upon parallel transport /irr'i simple bivectors are identified only up to spatial rotations. Only 
seems to remain manifest the coupling of the bivector norms, i.e the area of the triangles. 
This is the manifestation of 'ultralocality' in this geometrical setting, here due to the interplay 
between simplicity and parallel transport conditions induced by the non-commutativity of the 
star product. 

Note that the argument does not extend to general values of /9 in an obvious way, for the 
same reason that makes the constraint operator fail be a projector S^-^ S^- However it 
remains that upon commutation withl60|, a Lie algebra function gets conjugated not only by the 
holonomy /It-t'j but also by the Lagrangian multipliers of the simplicity functions. 

Arguing whether or not this feature is a serious problem from the point of view of quantum 
geometry is not our point here. Our point is to emphasize that it appears in our framework as 
an unavoidable feature following a clear geometrical construction of a dynamical theory for non- 
commutative tetrahedra. In fact, as it has been anticipated in the literature, the above argument 
shows that it is manifestly inherent to the non-commutativity of the bivector geometry, here 
entirely encoded into the star product. From the point of view of the canonical theory, where 
the tetrahedron states live in a tensor product of non-commutative spaces Ll{M.^), it is tied to 
the choice of quantization map. The presence of the star product thus encodes also quantum 
corrections, of which the above effect is a manifestation. 

This raises the question whether ultralocality survives in a semi-classical regime involving 
a commutative limit. The star product structure being dual to group composition, this limit 
corresponds to a linearization of the group. It can be formally defined by introducing a parameter 
e in the coordinates of the group manifold and to parametrize SU(2)^ group elements it = e*^" '^ 
for e.g by vectors Pu = \ smOn. The Fourier transform can be parametrized accordingly |36j : 
in the regime of small e, the S0(3) (resp. S0(4)) star product reduces to the usual pointwise 
product on functions of M"^ (resp. M^). The commutative regime should correspond by duality to 
the large spin limit of spin foam models. However since it amounts to linearize the holonomies, 
it could also be viewed as an analogue of the continuum limit in lattice gauge theories. Whether 
or not a proper development of the model around a commutative limit, at the level of the GFT 
|43j or its amplitudes, can be properly defined, and shown to tame the ultralocality feature, 
remains to be seen. 

5 Conclusions 

In this paper, we exploited a dual formulation of group field theories in terms of non-commutative 
bivector variables, which provide a duality between spin foam models and simplicial path inte- 
grals for constrained BF theories, to derive a new model for 4d gravity with Immirzi parame- 
ter. All geometrical variables remain explicit in this construction, which consists of inserting 
a constraint operator in a GFT for 4d S0(4) BF theory implementing the discrete simplicity 
constraints turning quantum simplicial BF theory into quantum simplicial gravity. 

Thanks to the framework chosen, we can keep the geometric content of all variables and of 
imposed constraints manifest at all stages of the construction. The resulting amplitudes for each 
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simplicial complex, generated in the Feynman expansion of the GFT, give a quantum simphcial 
version of the Holst-Plebanski formulation of gravity. We formulated these amplitudes both as 
BF simplicial path integrals with explicit non-commutative B variables and in terms of Wigner 
15j-symbols and fusion coefficients. The new model differs from existing ones in the literature; 
it imposes a different restriction on representation labels for quantum states. In particular, it 
does not lead any rationality condition for the Immirzi parameter. 

In light of this geometrical framework, we suggested a possible new perspective on the issue 
of the quantum correlations between neighboring simplices, often argued to be a problematic 
feature, for example, in the Barrett-Crane model. In our formalism, in fact, the relaxation of 
parallel transport condition is an unavoidable consequence of the very non-commutative nature 
of bivector variables and of the simplicity constraints, and tied to their quantization. Moreover, 
our framework is best suited for studying the geometric interpretation and consequences of this 
relaxation, as well as the semi-classical limit of the amplitudes. 

While this problematic issue certainly needs to be investigated further, we believe that the 
construction we performed suggests to consider the resulting model seriously as a candidate 
model for quantum gravity. In fact, its features appear all natural from the point of view of 
quantum geometry. Of course, this new model should now to be tested in all its aspects to 
support further or refute its validity. Note that its explicit formulation as a path integral for 
constrained BF theory, in contrast with the other existing models, should facilitate the study of 
its relation with a path integral quantization of continuum Holst-Plebanski gravity [31^ [32] . 

We also believe that the non-commutative formalism on which the construction is based, 
and that has proven useful already in different contexts and for different purposes, should itself 
be dwelled into in depth, to unravel even more aspects of (simplicial) quantum geometry. We 
have in mind the issue of symmetries, in particular the simplicial analogue of diffeomorphisms 
symmetry. This has been studied in the BF context in [H] and the analysis should now be 
extended to the 4d gravity model proposed here. Indeed, the bivector representation is the most 
suited one for defining such symmetries in a geometrically clear way. 

Finally, the really crucial question is whether it leads to an effective continuum dynamics 
of geometry, hopefully governed by some form of General Relativistic action, in the continuum 
limit. For this one has to study either the coarse graining of the lattice path integral appearing 
|44j in our amplitudes, or the renormalization flow and critical behaviour of our GFT model 
|41| . Also in this respect, our result offer a new, promising concrete model to analyze. 
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